Defects in Crystalline Packings of Twisted Filament Bundles: I. Continuum Theory of 

Disclinations 



O 
(N 

o 

Q 

00 



Gregory M. Grason 
Department of Polymer Science and Engineering, 
University of Massachusetts, Amherst, MA 01003, USA 
(Dated: December 30, 2011) 

We develop the theory of the coupling between in-plane order and out-of-plane geometry in 
twisted, two-dimensionally ordered filament bundles based on the non-linear continuum elasticity 
theory of columnar materials. We show that twisted textures of filament backbones necessarily 
introduce stresses into the cross-sectional packing of bundles and that these stresses are formally 
equivalent to the geometrically-induced stresses generated in thin elastic sheets that are forced to 
adopt spherical curvature. As in the case of crystalline order on curved membranes, geometrically- 
induced stresses couple elastically to the presence of topological defects in the in-plane order. We 
derive the effective theory of multiple disclination defects in the cross section of bundle with a fixed 
twist and show that above a critical degree of twist, one or more 5-fold disclinations is favored in the 
elastic energy ground state. We study the structure and energetics of multi-disclination packings 
based on models of equilibrium and non-equilibrium cross-sectional order. 
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I. INTRODUCTION 

Topological defects are crucial components of the 
structure and thermodynamics in many frustrated sys- 
tems in condensed matter [1!]. Perhaps the most well- 
known example of frustration-induced defects is exempli- 
fied by the Abrikosov phase of type-II superconductors 
where the application of a external magnetic field leads 
to the proliferation of vortices where the superconduct- 
ing order is destroyed and through which magnetic flux is 
threaded 0. One remarkable feature of such frustrated 
systems is appearance of defects in the ground state, in- 
dicating that the balance of competing thermodynamic 
forces requires highly heterogeneous distributions of en- 
ergy and structure. In geometrically frustrated systems, 
frustration arises intrinsically from the incompatibility 
of locally-preferred order with geometrical constraints on 
long-range ordering |3|] . Examples of geometrical frustra- 
tion abound in self-organized molecular and colloidal sys- 
tems. In bulk, three-dimensional materials, well known 
examples include liquid-crystal blue phases [3l as well as 
the twist-grain boundary phases of chiral smectics Q. 
When confined to two-dimensional surfaces, materials 
possessing a range of anisotropic order — nematic 
hexatic [?], smectic fT3, [HI, crystalline [ol [T2| - [l6j — are 
generically frustrated by the presence of non-zero Gaus- 
sian curvature, which couples favorably to the presence 
of disclinations in the ground states. 

Condensed phases of chiral polymers are subject to a 
unique frustration between two-dimensional (2D) order- 
ing perpendicular to the chain backbone and a micro- 
scopic tendency for intermolecular twist. Kamien and 
Nelson showed that in bulk columnar materials chiral- 
ity gives rise to two distinct types of grain-boundary 
phases — the tilt-grain boundary and moire phases — both 
composed of a rray s of dislocations lying in the plane of 
2D order [13, [l8| . Underlying the frustration of bulk 



columnar phases of chiral polymers is the geometrical 
coupling between in-plane displacements and filament 
tilts: inter-filament twist requires gradients in in-plane 
shear stresses along the molecular backbone. The frus- 
tration of chiral polymer assemblies is particularly impor- 
tant to the assembly of biological polymers, from extra- 
cellular proteins like collagen fiol and fibrin [l^ to the 
cytoskeletal filaments, f-actin |2ll. [23|. These are univer- 
sally helical molecules and are organized into densely- 
packed states in cells and tissues of living organisms [2j] . 
A number of studies show that the frustration between 
2D and chiral ordering gives rise to important thermody- 
namic properties including self-limiting assembly of bun- 
dles and fibers of chiral filaments |24| - [29| . 

Recently, we have demonstrated that frustration be- 
tween chiral order and 2D packing leads to another im- 
portant and surprising consequence, the restructuringof 
the ground state packing of twisted filament bundles [30| . 
In particular, we reported that twisted bundles are akin 
to 2D crystals formed on the spherically-curved sur- 
faces, in which out-of-plane geometry generates in-plane 
stresses that are partially screened in the ground state 
by topological defects, 5-fold disclinations, in the lattice 
order of bundle cross sections. Hence, it was proven that 
ground state in-plane order of sufficiently twisted bundles 
is irregular, possessing an excess of topological defects in 
the configurations that minimize the elastic energy of fil- 
ament packing. 

In this article we describe in detail the theory of 
disclinations in twisted filament bundles and its deriva- 
tion from non-linear continuum elasticity theory of 2D- 
ordered filament arrays. The unique symmetries exhib- 
ited by these materials lead to an intrinsic, non-linear 
coupling between in-plane strains and out-of-plane de- 
flections of the filament backbone: inter-filament twist 
makes it impossible to evenly space filaments in cross sec- 
tion. This analysis of the non-linear theory shows that 
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FIG. 1: The top (left) and side (right) view of a twisted 
filament bundle possessing and 5-fold disclination in the cen- 
ter. The filaments are colored to denote the tilt of filaments 
with respect to the center axes, with red filaments parallel 
the rotation axis and blue filaments most inclined. Here, the 
reduced twist corresponds to flR ~ 0.7 



the equations of mechanical equilibrium in bundles are 
subject a compatibility condition linking twist of filament 
backbones and in-plane stresses. In analogy to the Foppl- 
von Karman theory of elastic plates and membranes, this 
compatibility condition for filament bundles shows that 
twist generates in-plane stresses in bundles equivalent to 
those generated in a positively-curved elastic sheet pos- 
sessing a Gaussian curvature equal to Sfi^, where 27r/r2 
is the helical pitch of the bundle. Motivated by the fact 
that topological defects screen curvature- induced stresses 
in crystalline membranes, we derive the effective theory 
for disclinations in the cross section of twisted bundles. 
We show that above a critical of twist {flR)c = ■\/2/9, 5- 
fold disclinations are trapped in the cross section, where 
R is the bundle radius. An example of a twisted bun- 
dle ground state possessing an energetically-favorable 5- 
fold disclination is shown in Fig. [TJ At higher values of 
twist, multiple disclinations are stabilized in the cross 
section of crystalline bundle packings. Based on the con- 
tinuum theory, we study the discrete spectrum of elastic- 
energy minimizing configurations as a function of bun- 
dle twist which efficiently screen the buildup of twist- 
induced stresses. Further, we show that for a simple 
model surface-nucleated defects, the relatively weak re- 
pulsions between disclinations near the free boundaries 
of twisted bundles lead to a "pile up" of 5-fold defects 
just above the critical size for disclination stability. The 
excess defects entrapped within these non-equilibrium 
structures suggest that a limited mobility of chsclinations 
within the bulk of bundles drastically reduces the screen- 
ing of geometrically-induced stresses in "surface-grown" 
bundles. 

The remainder of this article is organized as follows. In 
Sec. |TT] describe the rotationally-invariant properties of 



the non-linear elasticity theory of 2D ordered filament 
arrays (i.e. columnar order). In Sec. IIIII we derive 
and study the equations of mechanical equilibrium for 
twisted bundles in the presence disclinations in the cross 
section, showing that one or more disclinations become 
stabilized by twisted- induced stress above a critical twist. 
In Sec. IIVI we present the results for multi-disclination 
cross sections in twisted bundles for both ground-state 
packings and a model of non-equilibrium bundle assem- 
bly. We conclude with a brief discussion of the surprising 
correspondence between crystalline order on spherically- 
curved surfaces and crystalline order of twisted bun- 
dles. Appendix A describes the coarse-graining of inter- 
filament interactions and the form of the non-linear elas- 
tic energy of multi-filament assemblies. An important 
result of this study is the exact and closed-form solution 
for the elastic energy of multiple disclinations in a 2D 
crystal, bounded by a free cylindrical surface. To our 
knowledge, no such solution has been previously pub- 
lished. As these results are more broadly applicable to 
finite-sized crystalline materials beyond the specific con- 
text of twisted filament bundles we include the full details 
of the derivation of the disclination-induced stresses and 
energies in Appendix B. In a subsequent paper we 
analyze the role of dislocations — "neutral", 5-7 pairs of 
disclinations — in ground-state order of twisted bundles 
as well as the structure and thermodynamics of multi- 
dislocation groundstates of twisted bundles. 



II. ROTATIONAL INVARIANCE AND THE 
NON-LINEAR STRAIN OF ORDERED 
FILAMENT ARRAYS 

In this section we construct the elastic energy that de- 
scribes deformations of an ordered array of filamentous 
elements. We assume the stress-free reference state of 
this energy to be an equidistant hexagonal array of uni- 
formly straight filaments, aligned along the z direction. 
To describe deformations, we consider the mapping of a 
material point filament array in its reference configura- 
tion, X — -|- zqz, to the position of the same point 
in the deformed configuration, r(x) = r^(x) -I- zz, where 
xj^ denotes the in-plane components x. We the define 
the in-plane displacement u(x) = r (x) — x^ , which the 
xy displacement of a point which, after deformation, is 
in a plane at z. Here, it is most convenient to adopt a 
"mixed representation" where we define u(x) as function 
of material coordinates, xj^, in the plane (a Lagrangian 
variable) and a function "current" vertical position, z (a 
Eulerian variable) |32| . 

As filamentous assemblies possess a type of two- 
dimensional order, it is essential to construct the strain 
in terms of a two-dimensional strain tensor that describes 
changes in inter-filament spacing in a plane perpendicular 
to the backbone orientation or chain tangent, t. The tan- 
gent is related to the displacement hy i ~ z + dzU± [ssj . 
This property is necessary to ensure that deformations 
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that shde filaments along their long axis cost no mechan- 
ical energy. The elastic energy of a material with this 
type of columnar order [s^ has the following form 

E=\ J dV{\ulk+2fiu,,u,,), (1) 

where Uij is the elastic strain tensor which has compo- 
nents in the x and y direction only. Here, A and fi are 
the Lame elastic coefficients characterizing the in-plane, 
isotropic elastic properties of the hexagonal filament lat- 
tice. The precise form of the strain tensor must be such 
that uniform rigid rotations of the reference state of the 
filament array require no elastic energy. The form of the 
strain satisfying these conditions is, 

Uij = ^ {diUj + djUt + diVL ■ djU ~ dzUidzUj) . (2) 

In addition to the standard symmetric derivatives in the 
linear elastic strain, u.y includes two non-linear contri- 
butions which ensure rotational invariance of the strain. 
The combination of the first 3 terms in eq. ([2]), famil- 
iar to the non-linear strain of 2D solids, is invariant 
under rotations around the z axis [s^]. The the final 
term is necessary to preserve rotational invariance un- 
der rotations around an axis in the xy plane [33|. It is 
straightforward to show, for example, that a rigid rota- 
tion around the y axis by an angle by 9 described by, 
Urj. = x{\ — secO) + ztanO leads to vanishing strain. 
This final non-linear contribution to the strain in eq. 
from the out-of-plane derivatives of u derives from 
the layered geometry of the filament array. The filament 
array can be decomposed into a series of rows, ruled sur- 
faces, stacked along directions in the plane perpendicular 
to the filament tangent. And not unlike the non- linear 
strain of fluid membrane stacks (i.e. smectics) [s^, the 
elastic energy is insensitive to sliding these rows along the 
tangent direction, and hence, the non-linear form of the 
strain ensures that the elastic energy measures distances 
perpendicular to t. That is, if neighbor filaments are 
separated by A in the xy plane, the non-linear strain is 
sensitive to the perpendicular projection of this distance, 
or 

Ai = A-t(t-A). (3) 

which measures the distance of closest approach between 
the filament pair. In Appendix A, for microscopic model 
in which filaments are treated as continuous curves whose 
segments interact via a central force law (such as a 
Lennard-Jones potential) we show that the potential en- 
ergy of the system is most naturally expressed as function 
of Aj^ alone. Furthermore, by explicitly constructing dis- 
tortions of this microscopic model from its mechanically- 
stable, hexagonally-ordered reference state, we find that 
the course-grained energy is identical to eq. ^ described 
by the non-linear elastic strain tensor, Uij, of the form eq. 

Thus, the contribution from —dzUidzUj ~ —titj to 
the form of the non-linear strain may be viewed as the 



result of the projection of inter-filament distance changes 
in a plane perpendicular to t, i .e. (5|Aj^p — 2uijSxi6xj. 

For the case of small strains, it is often sufficient to 
maintain the linearized strain tensor, Uij ~ {diUj + 
djUj)/2. However, as materials with columnar order are 
particular "soft" to tilt or bend deformations due to the 
lack of out-of-plane shear response, in-plane tilt deforma- 
tions become significant in many situations. It has been 
demonstrated that the anharmonic coupling between 
in-plane strains and filament tilt gives rise to out-of- 
planc buckling instabilities in columnar systems [33, 36]. 
In particular, when subject to a sufficiently large uni- 
form, uniaxial stress in the cross section of the material, 
the columns become unstable to certain long-wavelength 
modes along the long axes of the cylinder, which re- 
duce inter-filament spacing at the expense of out-of-plane 
bending. Hence, these filamentous and columnar mate- 
rials are subject to the same Helfrich-Hurault instability 
well known in smectic layer systems [40]. In the follow- 
ing, we show that this instability, is a particular example 
of a more general phenomena in which stress in the plane 
of filament order couples to out-of-plane deformations 
of flexible filament arrays. Due to the generic form of 
the non-linear coupling between filament tilt and elastic 
strain, we find that certain configurations of t necessarily 
introduce stress into the array, in a way that is formally 
quite similar to the relationship between in-plane stress 
and the Gaussian curvature of thin elastic sheets. For 
the purposes of the present study, we focus on the case 
of superhelically twisted bundles for which we demon- 
strate that the elastic energy ground state is, in general, 
highly complex and riddled with energy-minimizing con- 
figurations of topological defects in the lattice order of 
the bundle cross section. 

While it is important to keep in mind that is necessary 
to retain the contribution from diU ■ djU to the strain for 
large in-plane rotations that necessarily occur between 
distant planar cross sections of the twisted array, in the 
following we will drop the contribution from diU± ■ djU± 
when referring to Uij. Further, we use the fact that ~ 
dzUi to write the non-linear strain more simply as 

highlighting the geometrical coupling between in-plane 
strains and out-of-plane geometry unique to the non- 
linear theory of columnar materials. 

III. MECHANICAL EQUILIBRIUM IN 
TWISTED BUNDLES 

In this section we derive the equations governing the 
mechanical equilibria of filament bundles that are super- 
helically twisted. Our primary task is to demonstrate 
that a necessary consequence of this apparently simple 
geometry is a non-uniform distribution of the stress in 
the cross section of the bundle. Hence, the appearance of 
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stress in these bundles is an indication of the fact that it 
is geometrically impossibly to maintain a uniform spac- 
ing of filaments in the presence of twist [41]. Our for- 
mulation of the twist-induced stresses in bundles casts 
this geometric frustration in a surprising light. From the 
point of the view of the stress distribution, global twist 
of the bundle acts like uniform distribution of disclina- 
tions in the cross section whose density increases with 
the square of the bundle rotation rate. Thus, our second 
task in this section is to demonstrate the energetic cou- 
pling between the presence of twist-induced stress and 
topological defects in the lattice order of the filaments. 



A. Equations of mechanical equilibrium 

We consider a cylindrical bundle subject to a uniform 
rate, n, of twist around the central axis of the bundle. 
Additionally, we assume that the cross sections of the 
bundle are identical up to a rotation. In this case, we 
decompose the displacements into two steps: 1) u(x^), 
the in-plane displacement field from hexagonal order of 
a given reference plane of the bundle, say z = 0; and 
2) ua(xj_,z), the overall composite displacement result- 
ing from helical rotation of positions = xj_ -I- u at a 
constant rate, fl, around z along the bundle. 

uo(x) = [cos{Q.z) - l] {x'x + y'y) - sm{flz){y'x - x'y), 

(5) 

from which we find the in-plane projection of the filament 
tangent configuration. 



t^(r) ~ dzUn = f^r0. 



(6) 



This expression for t — given in terms the current po- 
sition, r, of the filament in cylindrical coordinates — 
describes a so-called "double-twist" texture, well-studied 
in the context of liquid-crystal blue phases f4!|. Here, fila- 
ment orientation rotates around the radial direction from 
the orientation, z, at the bundle center to a maximum 
tilt-angle. 



0max = arctan(ni?) ~ ^IR, 



(7) 



at the outer boundary, where R is the cylindrical radius 
of the bundle. 

In this paper, we consider the rate of bundle twist, 
il, to be quenched, say by the presence of strong, chi- 
ral inter-filament torques [27[ or, instead, by some ex- 
ternal mechanical force. According to helical symmetry 
described by eq. ([5]), the filament positions, displace- 
ments and strains are described by the state of a given 
xy plane at z = 0, with the bundle at other z described 
by uniform rotations of these configurations around the 
central axis. In this case, the problem reduces to solving 
for u(x) in the two-dimensional z — plane that min- 
imizes the energy per unit length described by eq. ([TJ, 
E/L = i J dA{Xu'^i, + 2fj,UijUij). Given the composite 



displacement field of eq. (O, the variation of the elastic 
energy with respect to components of u gives 



5{E/L) 
Sui 



where the stress tensor is simply 



(8) 



(9) 



The non-linear coupling between in-plane strain and out- 
of-plane filament tilt accounts for second contribution to 
the right-hand side of eq. ([S]). This represents a mechan- 
ical coupling of in-plane stress between different vertical 
layers of the bundle, which is absent from the linear the- 
ory of columnar order. 

Below we show that twist gives rise in-plane stresses, 
aij, or order (fiR)^. It is straightforward to show from 
eqs. ([5|) and ((6]) that the factor of in-plane tilt, ti and 
^-derivative of contribute to dz [tiO'ij] each an additional 
factor of (ili?)^ relative to the terms diatj in ([8]). Hence, 
for the case that small twist, we approximate the Euler- 
Lagrange equations, 6E/Sui = 0, by the order {flR)"^ 
terms and neglect the higher-order corrections deriv- 
ing from inter-plane coupling which contribute at order 
(rii?)'*. Based on this approximation, we arrive at the 
condition for mechanical equilibrium from the elasticity 
theory of isotropic materials, namely. 



diCTij 



0, 



(10) 



along with the condition of vanishing normal stress at 
the boundary. 



fiCTijir = R) = 0. 



(11) 



The problem of mechanical equilibrium in twisted bun- 
dles is one of 2D elasticity for which it is most convenient 
to solve aij in terms of an Airy stress function [37| . x, 
related to the stress by, 



(12) 



By construction, the stress derived from any function x 
is divergence- free, satisfying mechanical equilibrium ac- 
cording to eq. (jlOp . However, not every x corresponds 
to a physical configuration, which requires the existence 
of a displacement field u(x). That is, the solution for x 
must be compatible with the definition of strain, Uij. 

Similar to the well-known Foppl-von Karman theory of 
thin sheets iS^ , we derive the compatibility condition for 
the Airy stress from the constitutive relationship between 
stress and strain, 



ik^jidkdiX 
2jl 



AV2 



X 



(13) 



We cast the conditions that eq. (1131) can be solved to 
find a displacement field by taking the anti-symmetric 
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derivatives of strain, eik£jedkdeUij, yielding the following 
compatibility relation for the Airy stress 

(14) 

where we have used the definition of stain in eq. 
and i^o = 4^(A + a*)/(A + 2/i), is the 2D Young's mod- 
ulus of the array. For single- valued strain functions, the 
first term on the right-hand side must vanish, but in the 
presence of topological defects in the cross section, these 
terms generate non- vanishing sources for Airy stress. 

In general, two types of topological defects contribute 
to the Airy stress of the two-dimensional array: discli- 
nations, associated with singular configurations of 9q{'x.) 
the bond-angle field that points to the six-fold directions 
of lattice order; and dislocations, associated with singular 
configurations of u(x) that lead to no far-field rotation of 
the lattice directions The bond-angle field describes 
local rotations of the lattice directions in the cross sec- 
tion, and hence, is related to the anti-symmetric, in-plane 
derivatives of the displacement, 

6*6 = ]^e,jd^Uj. (15) 

In a hexagonal lattice, a disclination indicates a point 
around which increases or decreases by an integer mul- 
tiple of 27r/6, 

j) di-\7±0e = s, (16) 

where s = (27r/6)n is the topological charge of the discli- 
nation. From Stokes theorem we have the the areal den- 
sity of disclinations, s(x), 

s(x)=^s„<5(2)(xx-xi), (17) 

a 

where the sum over a indicates a sum over multiple discli- 
nations in the cross section. Likewise, a dislocation is 
defined in terms of a closed loop integral around which u 
changes by an integer multiple of the lattice spacing, a, 
along one of the six- fold directions, 

(j) di ■ \7 ^u, ^ b„ (18) 

where b is the Burger's vector. The areal density of dis- 
locations, b(x), is therefore, 

b(x)=^br<5(^)(x^-xl). (19) 

P 

By manipulating the displacement derivatives in the com- 
patibility relation it is straightforward to show ^15] that 
disclinations and dislocations generate point-like and 
dipole-like sources, respectively for Airy stress, 

7:<^ikej£dkde{diUj + djUi) = s(x) - V_l x b(x). (20) 



For the purposes of the present study, we focus on the 
case of disclinations only, b = 0, and solve exactly for 
the elastic energy of twisted bundles in the presence of 
arbitrary disclination configurations. In a subsequent 
study ^31j] , we take advantage of the fact that more com- 
plex topological defects, dislocations and grain bound- 
aries, may be constructed from neutral configurations of 
multiple disclinations, whose elastic stresses, according 
to the compatibility relation, are superposable. 

In addition to the defect-induced sources for Airy 
stress, eq. shows that certain gradients of in-plane 
filament tilt generate a more homogenous source for x 
due to the non-linear coupling of tilt and strain. We de- 
note these terms as the intrinsic twist, Kt, of the bundle, 
which can written as, 

Kt = -(-ike-jidkdptitj 

= X X ti)ti - (tx X Vi)ti] . (21) 

Indeed, by combining the disclination- and tilt-induced 
sources for Airy stress we arrive at the final form of com- 
patibility relation for x, 

-^Vlx = s(x) - Kt. (22) 

When Kt 7^ somewhere in the cross section of the 
bundle, eq. (|22p tells us that there must be in-plane 
stress. 

Though Kt is an unfamiliar geometric quantity, writ- 
ten in the form of eq. ([2T|) . we see that it is a total 
derivative of gradients of t in the xy plane reminiscent of 
the so-called "saddle-splay" of a vector field, n^. 

Kg = ■ [(Vi • nx)n± - (n^ ■ V i_)h4 , (23) 

a total derivative term in the Frank elastic energy of 
nematically-ordered materials [111 . The saddle-splay op- 
erator is most familiar when n is the normal to a surface 
that is weakly-defiected from the xy plane, say a mem- 
brane, in which case. Kg is the Gaussian curvature of 
the surface p^ . From this point of view, the intrinsic 
twist, Kt, plays precisely the role played by Gaussian 
curvature in the non-linear theory of thin elastic mem- 
branes. In membranes, the in-plane strain, couples 
to out-of-plane deflections [l^l, described by the height 
function, /i(x), according to 

~ i {d,uj + dju, + drhd.h) , (24) 

where the (m) in the subscript is used to distinguish the 
non-linear strain of membranes from the non-linear strain 
of filament arrays, eq. (|4]). The compatibility relation 
for the in-plane stress of membranes is known as the one 
of two Foppl-von Karman relations, describing the me- 
chanical equilibrium of thin plates [st']. This relation 
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was generalized by Seung and Nelson [l^ to include the 
presence of topological defects in the crystalline order of 
membrane and is identical to eq. (|22|) , with Kt replaced 
by the Gaussian curvature of the sheet, calculated from 
eq. (|23)) and fi_L ~ —\7±h. 

The compatibility equation for in-plane stress encodes 
two fundamental properties unique to the non-linear 
theory of membranes. First, it quantifies the in-plane 
stretching needed to deform intrinsically flat sheets into 
shapes corresponding to Kq ^ 0, in accordance with the 
renowned Theorema Egregium of Gauss [423. And sec- 
ond, it demonstrates the formal equivalence between the 
far-field in-plane stress distribution generated by local- 
ized regions of non-zero Gaussian curvature, and point- 
like, orientational defects in the crystalline order of mem- 
brane, 5- or 7-fold disclinations. Due to this second 
property, unlike crystalline order on planar 2D surfaces, 
topological defects — disclinations in particular — are well 
known to be necessary in the groundstates of crystalline 
order on surfaces with non-zero curvature, as disclina- 
tions of the appropriate charge are able to screen the 
curvature-induced stresses in-plane. An example of this 
frustration, important in structural studies of such di- 
verse materials as fullerenes [i^ , viruses [il] and particle- 
stabilized emulsions [l^, is known alternately as the 
Thomson problem, which seeks to describe the lowest- 
energy arrangement of repulsive, point-like particles on 
the surface of a sphere [45 1. 

Though ordered filament bundles possess a unique ge- 
ometry among 2D ordered materials, the mechanical frus- 
tration between in-plane stresses and out-of-plane geom- 
etry in these materials is remarkably similar to that pre- 
dicted by the non-linear elastic theory of membranes. In 
the case of filament bundles, the intrinsic twist plays the 
role play by Gaussian curvature in membranes, necessar- 
ily introducing distortions of the inter-filament spacing. 
Given the striking similarity, it is tempting, perhaps, to 
identify Kt with the Gaussian curvature of some surface 
or set of surfaces, say, to which filaments are normal. 
Such an identification easily fails for any twisted config- 
uration of filaments V_l x t_L 7^ 0, for which no such 
surfaces exist |46|, and yet Kt is demonstrably non-zero. 
For the case of super-helical twist considered here, we 
find, 

- ^{dltl + dltl - 2d,dytAy) = (25) 

Thus, superhelical double- twist generates a homogeneous 
and negative source of Airy stress in the cross section. In 
this sense, we find that the in-plane stresses induced by 
bundle twist are formally equivalent to those induced by 
positive Gaussian curvature in thin sheets. Formally, we 
may think of the intrinsic twist, Kt, as a uniform den- 
sity of negatively charged disclinations, continuously dis- 
tributed throughout the bundle cross section. Hence, as 
the bundle grows larger in radius, the integrated topolog- 
ical charge of this distribution grows as (fii?)^ ~ ^'maxj 
demonstrating that the frustration between uniform fila- 



ment spacing and bundle twist grows rapidly with bundle 
size. 

Finally, though twist generates stresses that are for- 
mally equivalent to those in spherically-curved mem- 
branes, it is important to point out a key difference be- 
tween the optimal packing of twisted bundles and the 
Thomson problem. Due to their closed topology, the 
bond network of point packings covering the sphere are 
required to possess exactly 12 excess 5-fold disclinations 
in all configurations, a requirement of topology rather 
than a ground-state property [l^ ■ Due to the presence of 
a free boundary, the problem of a twisted, 2D-ordered fil- 
ament bundle is more closely aligned with the problem of 
a finite-sized crystalline domain on a spherically-curved 
surface [J^ . The presence of a free boundary in of twisted 
bundle or curved crystal allows the net number of discli- 
nations (5- or 7-fold) to adjust according the energetic 
balance of disclinations needed to screen geometrically- 
induced stresses. Furthermore, as we show below, the 
requirement that normal components of stress vanish at 
the boundary of the domain leads to a strong-position 
dependence of the energetic costs for disclinations in var- 
ious positions in the bundle cross section. 



B. Elastic energy of disclinations in twisted 
bundles 



In this section we present the exact solution for the 
Airy stress of twisted, cylindrical bundles with multiple 
disclinations in the cross section, and derive the effective 
elastic energy depending only on fl and the positions and 
topological charges of disclinations. We seek a solution 
to the compatibility relation. 



1 

K'r 



Kj 



(26) 



The elastic energies of sufficiently twisted bundles show 
a clear preference for an excess of 5-fold disclinations, in- 
dicating that the groundstate packing of twisted bundles 
is complex and irregular. 

This basic result is clearly demonstrated by consid- 
ering the stress of a twisted bundle possessing a single, 
centered disclination. Denoting xt as the stress induced 
by intrinsic twist. 



Ko'^lxT = -Kt 



(27) 



we search for the form xt which satisfies stress-free 
boundary conditions. In polar coordinates the compo- 
nents of stress are related to x by. 



.-252 



'drX, 



and 



1700 — 9rXj 

<Jr4> = -dr{r~'^d^x)- 



(28) 
(29) 

(30) 
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Hence, for the case where xt is axisymmetric, the con- 
ditions that ar(ji = cTj-r = at the boundary reduce to 
R''^drXr=R = 0. The solution to eq. (f27| subject to this 
boundary condition is 

XT = ^(2i?V-rVT. (31) 

A single disclination at r = generates a stress denoted 
by Xsingie, Satisfying, 

i^o"'Vix..„<,/e -.sJ(2)(r), (32) 

and the stress-condition at r = i?. The solution for 
Xsingie for a Centered disclination in a 2D circular cross 
section is known from ref. |T2], 

Xs^ngle = (^) ^ [ l^l^/^) " 1/2]- (33) 

From the Airy stress functions we have the twist- and 
disclination-induced pressure fields, 

al, = ^{R'-2r')Kr, (34) 

and 

<r'^ = |7[l + 21n(r/i?)]. (35) 

Notice that the intrinsic twist induces as a dilation at the 
core of the bundle, cr^j. > for r < R/V^, and compres- 
sion at the periphery, a'^/^ < for r > i?/\/2. This trend 
is reversed for a positively charged disclination, such as 
a 5-fold defect, where cr^™^'^ < for r < Re~^/'^ and 
(7^™^''^ > for r > i?e~^/^. Hence, the induced stresses 
of 5-fold disclinations have a tendency to neutralize, or 
screen, the twist-induced stresses that derive from the 
non-linear coupling between in-plane strain and filament 
tih. 

As a result of this tendency, the elastic energy asso- 
ciated with bundle twist and disclinations are strongly 
coupled. From the Airy stress, the elastic energy can be 
computed from i J dVaijUij, which using eqs. (|9]) and 
P^ . can be manipulated to yield 

E/L=\j dA[K^\vlx? 

- ti~^eikejidkdi{dixdjx)Y (36) 

The second term in the integrand above is equivalent to 
the total derivative, fi~^dk{ukjdjx), whose contribution, 
therefore, vanishes due to the stress-free boundary condi- 
tion. Hence, the elastic energy derives directly from the 
in-plane pressure distribution, crfefe = V^^Xj 



The first term above represents the unscreened, non- 
linear elastic cost of bundle twist, whose fl- and R- 
dependence derives from the twist-induced strain, u^^ 
— (r2r)^/2. The second term represents the self-energy of 
a single disclination, which can be understood in terms 
of the constant axial strain required to remove or add 
angular wedges to the lattice, u^^ « s/{2tt). The fi- 
nal term in eq. (j37p represents the interaction between 
twist and disclinations. The sign of this term indicates 
a favorable interaction between twist and 5-fold discli- 
nations, s ~ -|-27r/6, due to mutually screening of the 
far-field stresses generated by these distortions. How- 
ever, while the elastic cost of twist depends continuously 
on n, disclination charge is discretized in units of 2tt/6. 
The minimal topological charge of 5-fold dislocations sets 
a critical value of bundle twist, 

(f7i?)e = y|~ 0.47, (38) 

above which the elastic screening by the disclination is 
sufficient to compensate for the defect self-energy cost. 
According continuum theory, bundles twisted beyond 
this threshold are unstable to the incorporation of one 
or more 5-fold disclinations in the cross section of the 
filament lattice. 

The calculation of the elastic energy of multi- 
disclination configurations follows a similar method to 
the case of a single, centered disclination. However, 
when one or more defects is off-center, the stress is no 
longer axisymmetric, requiring a somewhat more com- 
plex form. Though algebraically tedious, the solution to 
eq. (I22p proceeds by a standard multi-pole expansion of 
the disclination-induced stresses. Ultimately, the infinite 
series representing defect interactions with twist-induced 
stress as well as defect-defect interactions for disclina- 
tions in arbitrary positions can be resummed exactly to 
yield the relative compact expressions in the effective the- 
ory for defects in twisted bundles described below. As 
we know of no other published calculation of the exact 
disclination energy in cylindrical and circular crystals, we 
present the full details of this analysis in Appendix B. We 
note that this theory may have other applications, for 
example, for multi-disclination configurations in curved 
domains of 2D crystalline membranes. In the remainder 
of this section, we outline the solution and summarize 
the results. 

Consider a disclination at polar coordinates (p, (j) = 0) 
in the cross section of the bundle. Following ref. [4^ 
we decompose the Airy stress of the defect into a direct 
term, Xd, describing the singular stress at the core, and 
an induced term, Xi i due to the interaction of the defect 
with the free boundary. In the presence of a disclination 
of charge, s, the direct stress at a point (r, (j)) in the cross 
section has the form 

X, = ^r'Hnr', (39) 
where, r'^ = + — 2prcos0 and Kq^V^XcI = 
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s5^'^\r'). To determine the form of the induced stress, 
we consider the muhipole expansion of Xdi 



(40) 



Xd = {p^ ~ 2rpcos(f)) 

^ flnp-E,T=iH^)][cosM) r<p 
llnr-S;^^i i(^)"cos(n0) r>p 

The induced stress satisfies, V^Xi = for r < i? and has 
the generic form 



n=0 



[C„r-" cosM) + A.r-"+2 cos(n0) , (41) 



with Co = Ci = 0. The coefficients C„ and Z3„ are 
determined by boundary condition, = <^rr = at 
r = R. In terms of the Airy stress, the radial stress 
components vanish at the boundary when (48| 



and 



X\r=R = Kl COS0 + K2, 



R drX\r=R = Ki COS ( 



(42) 



(43) 



With the coefficients determined (see Appendix B), we 
evaluate the stress induced by the free boundary from. 



n=0 



4(n + l)Z?„r" cos(n(/ 



and from eq. (|B18|) we have the total pressure distribu- 
tion of the discHnation, 



'kk 



(^)- 



7^2 . 



Hf'/p)], (44) 
— 2rp cos (f) mea- 



where p = R^/p, and f'^ 
sures the square distance to an "image" defect of opposite 
charge outside of the bundle representing the screening 
of the elastic stress by the free boundary. 

As shown in Appendix B, the stress profile of multiple 
disclinations may be superposed along with the twist- 
induced stress, crjj,, and the total elastic energy may be 
calculated by area integration of (T^^(2i^o): yielding the 
effective energy quoted in ref. [soj , which may be decom- 
posed into 3 parts, 



E — Ett, 



E, 



disc 



(45) 



Here, Etwist/V — iKQ{VlR)'^ / 12^ is the non-linear elastic 
cost of bundle twist derived in eq. (1571) . From eqs. (IB34[) 
and (|B35|) we may combine the self-energy and twist- 
disclination interaction into a single disclination energy, 

Edisc7 



Ea^sclV = K,Y,vth^^Po.). 



(46) 



described in terms an effective, position-dependent po- 
tential for disclinations in twisted bundles. 



disc y^'-^Pa) — 227^ 



TT 



3(rji?)2 



2 s 2 



(47) 



Finally, the interactions between disclinations are de- 
scribed by 



(48) 



where 

^int (ra,r/3) 
|2 



167r2 



Pi 

i?2 



i?2 



■In 



{R^-pl){R^-pl)/R^ + \v^-vp\ 



(49) 



It is straightforward to check from eqs. (|T7)) and 
that both the elastic self-energy and the disclination in- 
teractions vanish when — > i?, a signature of the stress 
free boundary condition. 

The key result of this analysis is the mutual screen- 
ing of the elastic stresses of twist and disclinations in 
the bundle cross section described by the effective, single 
disclination potential V^l^^ {Q., pa)- Both the self-energy 
disclinations and the elastic interactions between twist 
and disclinations — respectively, the terms proportional 
to and —Sa{^RY in eq. (|T7)) — are strongest for discli- 
nations at the core of the bundle and decrease to zero as 
disclinations are pushed to the boundary of the bundle. 
However, while the self-energy always prefers to expel 
disclinations to the boundary, for the interaction between 
disclinations and twist can either be attractive or repul- 
sive, depending on the sign of Sa- As discussed above, 5- 
fold disclinations, with positive topological charge (sq — 
+27r/6), are attracted by the twist-induced stresses to 
the core of the bundle. For small twist, VI < fie, the self- 
energy of defects dominates over the elastic screening of 
twist, and vj^^J{n, Pa) > is repulsive, with a maximum 
for disclinations at the core, pa = 0. However, for suf- 
ficiently twisted bundles, 17 > f2c, the net elastic energy 
of single disclinations is attractive, favoring 5-fold discli- 
nations at the core the bundle. This attractive regime 
indicates that the screening of twist-induced stress dom- 
inates the elastic self-energy of defects, and ultimately, 
suggests that the optimal, ground-state packing of fila- 
ment bundles of fixed twist are irregular, including one 
or more topological defects in the cross section. 

Before analyzing the structure of these multi- 
disclination groundstates in the next section, we briefly 
discuss the global stability of 5-fold disclinations in bun- 
dles that are free to relax both the radial position of the 
defect, as well as the bundle twist. Fig. [5] shows a con- 
tour plot of the elastic energy, E{p, — Etwist + Edisc, 
of a bundle possessing a single 5-fold disclination as a 
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FIG. 2: The elastic energy landscape for a filament bundle 
possessing a single, 5-fold disclination in terms of defect ra- 
dius, p, and twist of bundle, SI. The color scale of the contour 
plot show lowest energies in purple and the highest energy 
values in red. Two locally stable minima are labeled: the 
global minimum, untwisted and defect free {p = R,Q. = Q), 
and the meta-stable twisted state with trapped defect at core 
{p = — {y/iR)~^). The dashed line shows a separatrix 
the divides the landscape into regions that flow along lines of 
force into each of the 2 minima. Defects originating to the 
right of this boundary are expelled, while defects originating 
to the left are pulled towards the core of the twisted bundle. 



function of radial position, p, and bundle twist, (fii?). 
In the absence of an external source for twist (i.e. im- 
posed torque) or intrinsic source for twist (i.e. chi- 
ral filament interactions), the global minimum occurs 
when disclinations are expelled from untwisted bundles, 
E{p = R^fl = 0) =0. Fig. [2] shows three additional 
critical points for < p < i? and flR > 0, where both 
the radial force on the disclination, fr ~ —drE, and the 
twist-moment on the bundle, M = —dnE, vanish. The 
point p = fl = is a local maximum describing the self- 
energy of a centered-disclination in an untwisted bundle, 
Eo = E{0,0) = KoV/2%%. Despite the absence an ex- 
phcit energetic preference for bundle twist, we find that 
there exists a meta-stable twisted state for p = and 
= ("v/Si?)^^, with a 5-fold disclination trapped at 
the bundle center for which £J(0,J7m) ~ Eq/A, the elas- 
tic cost of incorporating a disclination is significantly re- 
duced due to the screening disclination stress by twist. 
Finally, there is an unstable saddle-point separating the 
2 minima at = p = (1 — \/6/3y^^R and fig — flc, for 
which E{ps,fls) ^ Eo/3. 

The lateral motion of a disclination in 2D crystals re- 
quires a global repositioning of all lattice sites in the cross 
section, and therefore, we should expect the relaxation ki- 
netics associated with defects under the forces described 
by Vdisc to be extremely slow. Here we show that even in 



the absence of kinetic limitations on disclination motion, 
the non-linear elastic energy of filament bundles provides 
a thermodynamic mechanism to trap 5-fold disclinations 
in the cross section. Assuming a crude model of over- 
damped kinetics of bundle twist, dt^ — T^M, and de- 
fect motion, dtp = Ypfr, where Fji and Fp are damping 
coefficients, we expect weakly twisted bundles to expel 
disclinations that originate outside of the bundle core as 
the system evolves to a lower energy. However, when 
disclinations originate deep within the interior, the state 
of the bundle develops a finite twist, fi™ = (-y/Si?)"-^, 
and the disclination is entrapped within the bundle cen- 
ter. The boundary separating these two dynamical be- 
haviors is shown as dashed line in Fig. [21 While the 
entrapped disclination is ultimately thermodynamically 
unstable to expulsion from an untwisted bundle, there is 
an appreciable activation energy, i?o/12, proportional to 
the bundle volume, which must be overcome in order to 
remove the defect from the core of the twisted bundle. 
For this reason, we argue that despite the extensive elas- 
tic energy needed to form them, due to the elastic screen- 
ing provided by twist 5-fold disclinations formed at the 
core of two-dimensionally ordered filament bundles are 
especially long-lived, if not, permanent structural defects 
even in achiral materials, for which their is no intrinsic 
preference for bundle twist. 



IV. MULTI-DISCLINATION CROSS SECTIONS 
OF TWISTED BUNDLES 

In this section we analyze the elastic-energy ground 
states of twisted-filament bundles based on the contin- 
uum theory of topological defects derived in Appendix 
B and summarized in eqs. (|45l) - (H^ . These equations 
account only for the long-range elastic costs of deform- 
ing the array, and do not account, specifically, for the 
short-length scale effects captured by the core energy of 
disclinations, Ecore- In the following, however, we ne- 
glect this finite energy contributions in comparison to 
the elastic self-energy of disclinations that scales with 
Eseif ^ , and hence dominates in the limiting case 
of large bundles. 

A. Equilibrium packings of twisted bundles 

In this section, we consider the minimal- energy con- 
figurations constructed from an integer number of 5- 
fold, s = -|-27r/6, disclinations. According to the single- 
disclination elastic energy, eq. (|T7)) . for {Q.R) > {^lR)c = 
■y/2/9 5-fold disclinations are attracted to the core of a 
twisted bundle. While the non-linear stresses generated 
by twist attract disclinations to the bundle core, discli- 
nations of equal sign repel one another at long range, 
countering the tendency to entrap an arbitrarily large 
number of defects in the core. 

Determining the ground state of disclinations in 



10 



twisted bundles is, thus, not unlike the classical ground- 
state of long-range repulsive par ticles confined to finite, 
two-dimensional clusters |49l - l5l| . Such problems arise 
in a range of physical contexts including the Wigner- 
crystal states of electrons in mesoscopic particles [53j . 
2D-confined colloidal crystals [13] and multi- vortex states 
of superconducting discs [53| and quantum dots [SS^ . The 
models previously studied differ from the present specifi- 
cally in terms of the nature of both repulsive interactions 
and the effective confining potential. Unlike the very 
long-range elastic interactions between disclinations de- 
scribed by eq. (H^ . the long-range potentials studied in 
previous cases are typically inverse power law for or log- 
arithmic functions of separation. Additionally, the elas- 
tic screening by the stress-free boundary in the present 
model leads to disclination interactions, eq. (PU)) . that are 
weakened as either disclination approaches the boundary. 

The ground-state configurations of twisted bundles 
with only 5-fold disclinations were determined for twists 
in the range of < \ilR\ < 0.457r by direct comparison 
of the energies of defect configurations consisting of one 
or two disclination "shells" , a geometry common to few- 
particle groundstates of repulsive 2D clusters [i^. Each 
shell is composed of a C„ symmetric arrangement of n 
disclinations. For a given twist, the elastic energy of the 
bundle cross section is minimized over the radii of the in- 
ner and outer shells, as well as the angular correlation be- 
tween the shells. The energy, structure and total disclina- 
tion number of ground-state configurations are shown in 
FigM For bundle twist in the range ^/2/9 < < 1.05, 
the ground states possess a single shell of 1 < n < 5 discli- 
nations, while for higher twists, a second radial shell of 
disclinations is favored. 

While twist and its associated in-plane stresses in- 
creases continuously, the number of disclinations in the 
groundstate is quantized. Figure [3}3 shows the step- wise 
increase in disclination number in increasingly twisted 
bundles, a roughly linear function of (flR)^. The slope 
of the increase in disclination number is consistent, ap- 
proximately, with the number of 5-fold disclinations re- 
quired to "neutralize" the effective topological charge of 
defect and twist-induced stresses. From eq. (1^ . we see 
that out-of-plane twist leads to an effective topological 
disclination charge per unit area, —Kt = For a 

bundle of radius, i?, the number of s = -f27r/6 disclina- 
tions needed to neutralize this background charge can be 
estimated as. 



na{nR) 



KtA 
27r/6 



Q{nRf 



(50) 



which is plotted as dashed line in Fig. for compari- 
son to optimal disclination structures. While the slope 
of this curve is roughly consistent with the defect num- 
ber found by minimization of multi-disclination energy, 
notice that the number of disclinations predicted in the 
ground state is consistently less than uqIHR), suggesting 
that the neutrality condition is not precisely maintained 
for small bundles where much of twist-induced stresses 
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FIG. 3: In (a), we plot the elastic energy of ground state 
packings of twisted bundles predicted by continuum theory 
as a function of reduced twist. The distinct sections of the 
curve correspond to ground states with a distinct number of 
disclinations in cross section. Characteristic cross sections are 
shown for each disclination number, with points indicating 
the position of 5-fold disclinations in bundle. The dashed 
horizontal line in (a) denotes the upper bound on the elastic 
energy provided by an infinite Wigner lattice of disclinations. 
In (b), we plot the number of disclinations in the ground 
state in terms of bundle twist. Here, the straight dashed line 
corresponds to the number of disclinations needed for perfect 
"neutrality" of twist-induced stresses. 



are screened by the stress-free boundary condition. 

For large bundles, we expect the ground state predicted 
by the continuum theory of eq. (??) to converge asymp- 
totically a neutral, "Wigner lattice" of disclinations. The 
elastic energy of this configuration of the Wigner lattice 
can by calculated from the periodic solution to eq. ([^^ 
for the in-plane pressure for an infinite plane. 



(WL) 
'kk 



Kt 
i^o 



G^^^O 



iG : 



IGI 



(51) 



where G are the reciprocal lattice vectors of a 2D lattice 
of areal density Kx/ {2n/6) [5^. Assuming the ground 
state of the Wigner lattice is a hexagonal array of discli- 
nations, we calculate the elastic energy in the flR — >■ oo 
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limit from the lattice sum, 

i?(WL)/^ = ^ E ^ 0.00203X„-i. (52) 

Here, G = (27r/fi)[(2/i + fc)i/V3 + ky] where d = 

^J2tt /Z'^/'^Kt, the Wigner-lattice spacing of disclina- 
tions. This constant value is shown as a dashed, hori- 
zontal line in Fig. As indicated by the series of multi- 
disclination groundstates found for \{IR\ < 0.457r, the 
appearance of 5-fold disclinations in the cross section of 
twisted bundles accounts for a crossover from a (flR)'^ 
increase of the energy density to the ultimately twist- 
independent value of the infinite-radius, Wigner lattice 
state. 

It must be noted that in the forgoing analysis and 
derivation of the effect disclination energies, in the equa- 
tions of mechanical equilibrium, eq. (jlOp. we have ne- 
glected terms that are higher-order in (fii?)^ that derive 
from coupling of stresses in different vertical layers of 
the bundle. Clearly, the inclusion of these effects must 
have a significant impact on the structure and energy 
of multi-disclination ground states for large values of re- 
duced twist, when {IR\ 1. However, the predicted 
value of {flR)c — 0.47 where the bundle becomes un- 
stable to a single disclination, as well as the transition 
to multi-disclination ground states, fall well within the 
regime of twist where we expect the corrections to the 
approximate theory to be modest. 



B. Non-equilibrium disclination arrays in 
surface-assembled bundles 

The motion of disclination lines in twisted bundles of 
filaments is slow, and in many realistic scenarios the po- 
sitions of disclinations that form in the cross section of 
twisted bundles are likely to be quenched during the as- 
sembly process. The displacement of a disclination, say 
by a distance of one lattice spacing in the cross sec- 
tion, requires the absorption or emission of a dislocation, 
a neutral disclination pair, which in turn must diffuse 
from the disclination to the free boundary of the crys- 
tal [53, [Ell. Hence, the motion of a disclination over 
a distance of order R, requires a sequence of roughly 
R/a dislocations to glide a distance of order R/a lat- 
tice spacings. We may crudely estimate the timescale for 
disclination motions within a 2D crystal grow at least 
as quickly {R/a)'^L, where the factor of bundle length L 
accounts for the fact that disclinations run along the full 
length the bundle. All of this suggests that over the times 
scales during which the self-assembly process occurs, it 
is quite reasonable to expect that disclination motion is 
negligible, and therefore, the most probable state of even 
highly-twisted filaments assemblies will be markedly dif- 
ferent from the equilibrium, multi-disclination states dis- 
cussed in the previous section. 



In this section we explore a simple model for the 
non-equilibrium formation of twisted-filament bundles in 
which the topological reorganization of the cross section 
is inhibited by the sluggish dynamics of disclinations. In 
this model we assume that disclinations can only form at 
the free surface of a bundle as new filaments adhere to 
a partially-formed bundle. Upon nucleation of disclina- 
tions within a surface layer of newly-adsorbed filaments, 
defects are then assumed to be frozen in position as fur- 
ther layers of absorbed filaments are assembled on the 
bundle. 

To study the cross sectional order within this model 
of surface assembly we construct the following algorithm 
where concentric radial layers of filaments are absorbed 
sequentially to a bundle of fixed rotation rate, fl. 

1. Add radial filament layer of width a to pre-existing 
bundle, increasing outer filament tilt by an incre- 
ment, 69 = afl. 

2. Consider addition of new disclinations at radial 
position, R — a, by minimizing eqs. (j45|) over 
new defect number and angular position. Addi- 
tional disclinations are subject to repulsive interac- 
tions within incipient layer, as well as with defects 
trapped within bundle from earlier stage of growth. 

3. Energetically favorable defects are trapped within 
growing bundle and maintain their in-plane posi- 
tion upon growth of further surface layers. 

The minimization over disclination number and angular 
position is carried out numerically, and at each step in 
the radial-assembly process, the total number of disclina- 
tions, ndisc, and total energy of the twisted-bundle cross 
section are determined. 

In Figure 2] we show results for the surface- assembly 
model of twisted bundles with 5-fold disclinations in the 
cross section. These demonstrate that structures of the 
multi-disclination cross sections are both very distinct 
from the equilibrium cross sections of Fig. [31 as well as 
crucially sensitive to the tilt per filament layer, 66 = fla. 
According to eq. (|T7)) . the stability of disclinations in 
twisted bundles is determined entirely by the reduced 
twist, flR, the critical value of which is a/2/9 ~ 0.47, 
is the same for all inter-filament spacings a, allowing us 
to correlate the magnitude of the tilt increment with the 
size of a bundle at the critical twist. A smaller value of 
tilt per layer, Qa, implies that the bundle will be larger 
when it reaches the critical twist for stable disclinations. 
This critical radius for disclinations is R^,/a = 0.47/(56', 
in units of filament spacing. 

Figure shows the energy density of surface- 
assembled bundles for four values of tilt increment: 
60 = 0.01,0.02,0.05 and 0.07. Surface-assembled bun- 
dles share the same (fii?)^ energy-dependence below the 
critical value of twist, and each show a characteristic 
maximum in E/V due to the presence of 5-fold discli- 
nations at higher twist. Notably, these energies "over- 
shoot" the ground-state energy density above the critical 
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FIG. 4: The numerical results of the non-equilibrium model for radial growth of twisted bundles. In (a), we show the elastic 
energy vs. twist for bundles grown with for different values of tilt increment: 59 — 0.01 (blue), 0.02 (red), 0.05 (orange) and 
0.07 (green), (b) shows the total number of disclinations in the bundle cross-section for these same values of SO, with cartoons 
depicted the position of 5- fold disclinations at the upper limit of fii?. In (c), the number of disclinations in the first radial 
"shell" (the first pile-up), is plotted as function of tilt increment. 



value of riR, and the extent of excess energy (relative 
to the groundstate) depends sensitively on tilt increment 
per layer, increasing dramatically as 69 is diminished. In 
this model of assembly, larger bundles, measured in terms 
of i?*/a, screen twist- induced stresses less effectively than 
smaller bundles by the incorporation of 5-fold defects at 
their free boundaries. This effect can be traced to the 
number of disclinations which are captured in the first 
few radial layers of growth beyond the critical size. 

The number of disclinations captured, ndisc, at differ- 
ent values of VlR, is shown in Figure. for the same 
four values of 69. Each assembly exhibits an immedi- 
ate "pile-up" of 5-fold disclinations in the outer layers as 
soon as the bundle grows beyond the size necessary to 
stabilize a single disclination, R,, — y^2/9Q^^. That is, 
just above the critical size, these assemblies jump from 
a defect-free state to assemblies with n^isc = 11,6,5 and 
3 disclinations, respectively, for 69 = 0.01,0.02,0.05 and 
0.07 (osj. Hence, in comparison to energy-minimizing 
configurations of defects, when surface-assembled bun- 
dles reach the critical size, the assembly process tends to 
overcompensate for twist-induced stresses with a num- 
ber of disclinations that far exceeds the number needed 
to neutralize the twist-induced stresses. The number of 
excess 5-fold defects correlates directly with the excess 
elastic energy at large twist (Fig. |3]a) as the defects are 
eventually incorporated in the core of the bundle where 
the elastic costs of inter-defect repulsions are more costly. 
As shown in Fig. |3]b the disclination pile-up tends to lead 
to the formation of radial bands of concentrated disclina- 
tions. The larger the number of disclinations in the first 
pile-up, ndisc(l), the higher the value of reduced twist, 
rii?, needed to produce the secondary band, as the elas- 
tic repulsion with the initial shell tends to stabilize a 
growing bundle to defects at further radial shells. 

In Figure |4] (c) we analyze the number of disclinations 
in the first radial "shell", ndisc(l), as a function of 69. 
Here ndisc (1) is identified from the initial plateau value 
of ndisc beyond the critical twist. Two trends are imme- 



diately apparent. First, the value of ndisc (1) is oscillatory 
with respect to 69. Second, the mean number disclina- 
tions (the amplitude of the oscillatory envelope) in the 
first shell decreases with 69. The oscillations in ndisc (1) 
are simply due to the commensurability of a finite num- 
ber of layers, Ng, of finite radial thickness, a = 69/Vl, 
with the critical bundle size, i?*. The number of layers 
of radial growth achieved when the bundle has grown just 
beyond the critical size is simply iV* = 1 -I- int[ ^2/9/(56'] . 
In the model of radial growth, a bundle will therefore be 
quenched just beyond the critical twist by an amount, 
6{nR) = NJ9 - v/279, at the point when the first shell 
of disclinations forms. It can shown that due the integral 
nature of layer addition, this degree of overtwist oscillates 
with increasing 69 over a range and 69 with a period 
roughly corresponding to 69"^ / a/2/9. 

The decrease of the mean value ndisc (1) with 69 can 
be related to the weakening of repulsive interactions of 
disclinations near to the free boundary of the bundle. Ac- 
cording to the effective energy of disclinations in twisted 
bundles, eq. P5)) . ndisc (1) defects added near to the bun- 
dle boundary at p = R — a leads to an elastic energy gain 
of roughly, —n^isc{^)KQ5{^lR)a?. Opposing the incor- 
poration of multiple disclinations in the surface layer is 
the cost of inter-disclination repulsions. From eq. (|49|) . 
when defects approach the surface of the bundle such 
that Pa ^ R (or pfj — ^ R) the inter-defect potential ap- 
proaches zero as 



2 ^2 



lim Fj„t(rQ,r^) 



[R^-plfjR^^p l) 
327r2i?6|r„-r/3| 



(53) 



Assuming a multi-disclination shell with ndisc (1) defects 
spread evenly within a radial layer at p — R — a, we esti- 
mate the total cost of disclination repulsions per shell as 
roughly n'^^^^{l)KoV{a/ R)'^ . Although the favorable cou- 
pling between disclinations and twist vanishes as defects 
approach the bundle surface like (a/i?)^, the strength of 
disclination-disclination repulsion vanishes more rapidly 
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as (a/R)'^, and thus, provides a weak resistance to the 
pile-up of dischnations at the surface layer of the bundle. 
Optimizing the net elastic energy gain over the number 
of dischnations in the first shell we expect, 

where the overtwist may be approximated by S{nR) ~ 
66. This analysis shows, consistent with numerical re- 
sults, that the number of dischnations in the first pile-up 
layer diverges as (56' — 0. 

The results of the simple kinetic model of bun- 
dle growth suggest that the efficiency of twist screen- 
ing by elastically favorable 5-fold dischnations in non- 
equilibrium bundles is highly dependent on tilt per radial 
layer, 69 — fia. Interactions between surface defects in 
bundles composed of relatively few filaments upon reach- 
ing the critical twist (i.e. (i?*/a)^ = 69^ is relatively 
small) tend stabilize the cross section against pile-up of 
excess dischnations. While the excess defects trapped at 
the surface of relatively larger bundles ultimately lead 
to additional elastic costs to further bundle growth as 
the surface defects must ultimately be incorporated into 
the core of the bundle, where inter-disclination repulsive 
costs are greatest. 

V. DISCUSSION 

In this article, we studied the infiuence of out-of-plane 
geometry on the in-plane packing of twisted and two- 
dimensionally ordered filament bundles based on non- 
linear continuum elasticity theory. Formally, this influ- 
ence derives from the unique form of the non-linear 2D 
strain tensor describing deformation of columnar order, 
which couples inter-fllament strains and tilting of flla- 
ment backbones. Analysis of the equations of mechan- 
ical equilibrium in filament bundles shows that helical 
twist generates stresses in the bundle cross section that 
are formally equivalent to those generated by positive 
(spherical) Gaussian curvature of the 2D elastic sheet. 
In analogy to the problem of the curved, crystalline 
membranes, we also find here that these geometrically- 
induced stresses may be screened, in part, by 5-fold disch- 
nations defects in the cross section of the bundle. 

We derived the effective theory of dischnations in the 
cross sections of twisted bundles and, based on this the- 
ory, studied the multi-disclination states of bundles for 
fixed reduced twist, (QR). In equilibrium structures, 
dischnations appear above a critical twist, {flR)c = 
•\/2/9, and unlike the energy of the defect-free bun- 
dle, which grows unbounded with $7*, the energy of 
the disclination-possessing groundstates remains below 
a finite, upper bound provided by the infinite- twist, 
Wigner lattice of dischnations. In a simple model of non- 
equilibrium assembly of the twisted bundles, we show 
that the strong screening of disclination-disclination re- 
pulsions by the free boundary of the bundle leads to an 



instability to "pile-up" excess dischnations at the sur- 
face of a bundle when the bundle grows to the critical 
value of reduced twist. As the number of defects in these 
non-equilibrium packings far exceeds the number needed 
to neutralize the geometrically-induced stresses, in this 
model of bundle growth the presence of 5-fold dischna- 
tions is much less effective at reducing the in-plane elastic 
stresses that may act to limit the lateral growth of bun- 
dles during the assembly process. 

Underlying complex structure of the cross-sectional or- 
der in twisted bundles is the surprising and unusual con- 
nection highlighted by non-linear elasticity theory be- 
tween crystalline order on spherically-curved surfaces and 
in twisted bundles. We note here that the results of two 
similar models of chiral and biaxially-ordered materials 
anticipate this connection. The first was demonstrated 
by Kamien for a chiral, liquid-crystal phase with both ne- 
matic and hexatic order, the so-called N + 6 phase [soj . 
Here, the six-fold bond-orientational order is represented 
by a vector field, N, in a plane normal to the local ne- 
matic director, n. The geometric relationship between N 
and n — encoded in the Mermin-Ho relation — leads to the 
impossibility of uniform hexatic order for certain textures 
of the nematic director. Indeed, for the double-twisted 
texture studied here, the intrinsic geometry of the ne- 
matic frustrates in-plane hexatic order in precisely the 
same way that curvature frustrates hexatic order formed 
on a spherical surfce. Hence, in large double-twisted do- 
mains of the -|- 6 phase, just as in spherical hexatics, 
positively charged, 5-fold dischnations are favorable to 
screen geometrically-induced stresses. 

Further evidence of the formal equivalence between 
frustrated order in twisted bundles and on spherical sur- 
faces was demonstrated by Kleman in studies of ideal 
packings of twisted filaments in an intrinsically curved, 
three-dimensional space, namely [13, [Ml- The results 
of these studies derive from the Hopf fibration, which 
decomposes into linked geodesic curves, great cir- 
cles, that are each associated with a unique point on 
[31 . The geometry of has the property that: 1) 
any two great circles are equidistant and twist helically 
along their length; and 2) the distance between them in 
is half of the separation of their associated points on 
S^. Hence, any equally-spaced arrangement of points on 
S'^ corresponds to a perfectly-packed (maximum density) 
and twisted configuration of curves, or filaments, in S^. 
In this case, the problem of point packings on has a 
one-to-one mapping onto a twisted fiber packing of S'^ . It 
is less clear at present, however, how this mapping may 
be used to study twisted filament packings in Euclidean 
space, as no projection from to R^ preserves curve 
separations, leading to a second level of frustration as- 
sociated with the flattening out of ideal packings of 
in addition to the frustration of point packings on al- 
ready encoded in the flbration. As with a certain other 
geometrically frustrated problems — notably liquid crys- 
tal blue phase textures and icosahedral order in amor- 
phous materials 9|~exploring the structure and energet- 
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ICS of ideal, twisted filament packings of projected to 
may ultimately provide a useful alternative avenue 
for studying chiral filament packing from a fundamen- 
tally geometrical point of view f62| . 

We conclude this article with a few remarks regard- 
ing the potential for observations of multi-disclination 
ground states in filamentous systems known to form 
twisted bundles or fibers. According to continuum the- 
ory, in both equilibrium and non-equilibrium bundles 
the key parameter determining the appearance of 5- 
fold disclinations is the reduced twist, flR, with a crit- 
ical value (r2i?)c — \/2/9 ~ 0.47. This corresponds 
to a maximum tilt angle of the outer most filament of 
6c — arctan[(f2_R)c] — 25.2° with respect to the bundle 
axis. This high degree of bundle twist has been demon- 
strated in coarse-grained simulation models of filament 
assembly [G^, and hence 5-fold defects are readily ob- 
served in the cross sections of simulated assemblies of 
chiral filaments [1^. However, this critical tilt exceeds 
the value observed for all if not most twisted, biofilament 
assemblies. The collagen fibers that are observed to twist 
helically reach a maximum tilt angle of 15 — 16° [l^. The 
pitch and radius distributions of fibrin fibers have been 
well characterized in vitro by Weisel et al. (23 |. showing 
a maximum radius of roughly 50 nm and a mean pitch 
1.9 /im, corresponding to an outer filament tilt of ap- 
proximately 9° , falling well below the threshold at which 
disclinations become stable. It is notable, however, that 
under certain conditions fibrin fibers are observed to grow 
to much larger radii, at which point the in-plane pack- 
ing becomes highly-disrupted, leading to a significantly 
more diffuse, and possibly fractal, organization at the 
fiber periphery. In light of the present study, we specu- 
late that the appearance of disrupted in-plane packing of 
fibrin fibers beyond a critical size may be a remnant of 
the disclination pile-up that is predicted to occur in the 
non-equilibrium assembly model. 
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Appendix A: Coarse-Grained Elastic Energy of a 
Microscopic Model of Self- Assembled Filament 
Bundles 

In this appendix we analyze a microscopic model of fil- 
amentous assemblies and deduce the form of the elastic 
energy via a coarse-graining of the microscopic degrees of 
freedom. Our purpose is to demonstrate explicitly that 
the form of the elastic energy in eq. ([1]) described specif- 
ically in terms of the non-linear strain of eq. ^ arises 



naturally for ordered filament bundles forming under the 
infiuence of generic and mutually-adhesive interactions. 
For simplicity, we assume that interactions between a 
pair of filaments, labeled a and /?, may be described by 
the following interaction energy, 

Uaf3^^ J dSa j dsp v(\Ya{Sa) - Vp{sp)\), (Al) 

where J dsa (or J dsp) denotes the integration over the 
arc-length of the central backbone of filament a (or /3) 
which is described by ra{sa) (or r^(s^)). Here, v{r) is 
the pair potential describing interactions between seg- 
ments on different filaments, which we assume to be 
short-ranged and to be a function of for simplicity. We 
carry out the second integration over sp for a fixed Sa by 
identifying as the arc-position on /3 at the same ver- 
tical height as ra{sa)- We may then approximate r^(s^) 
to second order in Ss = sp — as. 

Map) = Ms'J + WJSs + ^ii^npis'JiSsr + ■■■, 

(A2) 

where t/3(s„), np{s'^) and Kpi^s'^^) are the tangent, normal 
and curvature of filament /3 evaluated at . At each po- 
sition Sa we denote the in-plane separation between the 
curves as Aa^(sa) — ra{sa) — r^(s^), while the square 
separation between two points at different vertical posi- 
tions is approximately given by, 

|Ar„^|2 = A^^-2A„/3-t;3 Ss 

+ {1-K0np-Aap) 6s'' + 0(6.8''), (A3) 

where we have suppressed the Sq dependence for clarity. 
In densely packed array, we expect filaments to be suf- 
ficiently straight on the scale of interfilament separation 
so that K/3n^ • Aq.^ <C 1 and we may ignore the curvature 
correction to |ArQ/3p above. With this approximation, 
the separation between curves a and (3 obtains a mini- 
mum for arc-separation, 

6s^^Aap-ifi, (A4) 

which denotes the closest point to ra{sa) on r/3(s/3), or 
the distance of closest approach, 

Alp^^ EE \Arap\l = A% - [Aap ■ ipf. (A5) 

Hence, Aq,0,_l is perpendicular to tp at sp = s'^ + <^s*, 
and expanding around the distance of closest approach 
As — Ss — Ssif, we have the separation between points on 

a and /3, 

\Arap\'^Alp,^ + {Asr. (A6) 

Inserting this into eq. (jAl[) and integrating over As we 
find that we may express the interaction between a and 
/? as the integral over Sq of a natural function of A^^ _|_, 
the local distance of closest approach, 

C/a/3 = / dSaVeff{Aap,±), (A7) 
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where taking the length of the filaments to be infinite 

V,ff{A^) = r du «(|A^|v/TT^). (A8) 

For example, for filaments whose segments interact via a 
Lennard- Jones potential, v^^'^\r) = e[(cr/r)^^ — (cr/r)^] 
it is straightforward to show that, 



7120" 
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76 CT 



(A9) 



where 7„ = /:^^dV(l + ^^'r/'- 

Provided that the interactions are sufHciently short- 
ranged, it is clear that pair-wise interactions between 
infinitely long filaments can always be decomposed, as 
in eq. (|A7p . into the arc- length integral along a given 
curve a of an effective potential energy. This effect po- 
tentail depends only on Aq,/3,_l = A^^ — ■ Aq/3, the 
local nearest distance between a and /3. Assuming cohe- 
sive interactions between filaments, ^^//(Aj^), obtains a 
minimum for some value, | A^| ~ Aq. Expanding around 
this separation we will have. 



(AlO) 



with T4//(Ao) — —e and k — (PVef f / d/\\. We assume 
that interactions are sufficiently short-ranged that we 
may truncate interactions beyond the nearest neighbors, 
and the ground state is a uniformly-spaced, hexagonal 
array of straight filaments. Deformations from this state 
may be described by the energy. 



E = 



dSa 

dz 



'a/3, 



Ao) 



(All) 



where Eq is the energy of the equally-space, reference 
configuration and the sum is carried out over nearest 
neighbor pairs in the hexagonal array. Here, dsa/dz = 
{ta-z)~^ due to the change in arc-length within a vertical 
layer of width, dz. For a given filament pair that is ini- 
tially separated by the nearest-neighbor vector AoTq^, we 
may write the separation in the distorted state in terms 
of two continuous functions, the displacement, u(x), and 
the filament tangent, t(x) ~ z + 9zU, 

Aa^^j^ ~ AofQ,3-)-Ao(fa/3-V)u-Ao(fa/J-l9zU)t, (A12) 

where we have retained corrections of first order in deriva- 
tives of u, the in-plane displacement. The change in the 
square-separation is straightforward to compute, 



|AJ2-A^ 



Agfifj (2diUj + diU ■ djU — dzUid, 



(A13) 

where we have suppressed the indices, a and on Tap 
for clarity. Substituting into eq. (jAll[) and using the 
fact that for a given lattice position, a, the average over 
nearest neighbor vectors in the lattice yields, 

9 

(nfjfkfi) = -^{SijSke + 5ik5ji + SuSjk), (A14) 



and we can write the change of energy upon deformation 



E~Eo^^J dV{Xulk + fiu,,u,, } . (A15) 



Here, the elastic constants are 



9kpo 



(A16) 



where po is the in-plane density of the reference configu- 
ration and the strain tensor. 



diU ■ djU— dzUidzUj) . (A17) 



Notice that Uij is "mixed representation" strain |32| . 
This means that the u is implicitly constructed as a func- 
tion of the undeformed, reference state material coordi- 
nates for in-plane variables x and y, a Langrangian de- 
scription, while it is a function deformed positions for the 
vertical z coordinate as in an Eulerian description (40j. 



Appendix B: Elastic energy of multi-disclination 
configurations by multipole expansion 

1. Single disclination stress and self-energy 

In this section we present the solution of the Airy stress 
for a single disclination at a radial position, p, in the cross 
section of the twisted cylindrical bundle of radius R by 
multipole expansion. The infinite series that result from 
these expansions may be resummed, making use of the 
following series for \z\ < 1, 



n=l 



In(l-z), 



(Bl) 



■a 

5: = -1 - z-Hn(l - z), (B2) 



E '-^^ = ln(l -z) + U,{z), (B3) 

n=l 



Z 

E = l + z-' HI -z)- ln(l - z), (B4) 

n— 1 ^ ^ 



V f 

^ ri^ln + 1) 

n— 1 ^ ' 

and 

oo 

E 



■nP'{n — 1) 



= -1 - ln(l -z)+ ln(l -z) + Li2(z), 

(B5) 

2z+ln(l-z)-zln(l-z)-Li2(z), (B6) 
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where Li2{z) = X^J^i = jv? is a polylogarithm. 

The coefficients of the induced stress, Xi^ are chosen 
to satisfy the stress free conditions, eqs. (|^ and (^5]). 
From eq. we have the coefficients of the multi-pole 
expansion of the direct stress, 



Xd 



where, 



1 E,T=o xi"> cos(n(/i) r > p 



Xd,< = -2prlnp- p2 _ _^ 



9 ? 



n{n — 1) n(n + 1) 



and 



n(n — 1) n(n + 1) 



(B7) 



(B8) 



(B9) 



(BIO) 



(Ell) 



(B12) 



(B13) 



For the case of stress-free boundaries, the elastic en- 
ergy only depends on (7^^ = V5_X; and hence, the terms 
in eq. (|4T|) proportional to the coefficients C„ do not con- 
tribute to the energy. Defining D„ = {KQs/8n)dn, from 
R drx\r=R = Ki coscf) we find that. 



do 



InR 



1 



2 2i?2' 

From x\r=R - R drx\r=B^ = K2 we find, 

^3 



di 



_P_ _ _P_ 

i?2 2i?4 ■ 



(B14) 



(B15) 



Finally, the as all terms proportional to cos(n(/)) in x and 
drX vanish at r = i?,, wc find the condition 



R drX^^l\r=R - nx';i\r=R + 2d„i?"+^ = 0, > 2) 



which is satisfied by 



i?"+2 \R 



i?2 



n rt + 1 



(B16) 
(B17) 



From these coefficients, we may compute the induced 
pressure, cr^^ = V5_X*, 



4fc = ^ XI "^(^ + 1)^"''" cos(n0) 



Stt 
Stt 

Kos ( 



n=0 



41ni?- 2 



R' + p^ 
i?2 



^ li?2'^^ y V i?2 



Stt 



-41ni?-4AT 



2 71-2 



i?2 



!ln(f'2/p2)|, (B18) 



where p = R^/p and f'^ = p^ + r"^ — 2rpcos(f>, and we 
have used eq. (jBl|) and the geometric series. Note the 
appearance of a singularity at f' = 0, which indicates 
the presence of an "image" defect outside of the bun- 
dle boundary, at a distance p > R from the center of 
the bundle, which screens the normal stresses at the free 
boundary. 

The next task is to use the disclination-induced stress 
to calculated the elastic self-energy of disclinations and 
the elastic interaction between defect- and twist-induced 
stresses. Due to the orthogonality of multipolc terms 
with respect to (j) integration, to perform area integration 
of the elastic energy it is most convenient to work with 
a full expansion of the defect induced pressure, crf^'"^ = 



disc 
"kk 



where 



Stt 



=0 ^ cos(n0) r<p 

(n) . / ,N . ' y^^y) 



Eoo (n) / j,\ ^ ' 

„=o cos(n0) r> p 



Un{p/R) 



'i?2' 



(B20) 



1 /?'\" {n + 1) f rp 







L n 





and 



4°)=4In(r/i?) + 2-2|j, 



(B22) 



(n>l) _ 



Defining 



1 (n -f- 1) / rp\" 

frp_Y/p_^^ 



\R^J \R 



4") = r dr r{a^^\r)y 
Jq 



(B23) 



(B24) 
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and 



Using these expressions and eqs. (|B1[) - (|B6[) we find, 



E 



in) 



(B25) 



we may write the elastic self-energy of the disclination as 



^ oo 



(") 
< 



i?2 



- 1 



8p2ln(i?/p). (B33) 



Eself /L 



647r2 



(0) , ^(0) 



- 

n=l 



(B26) 



Evaluation of the radial integrals is straightforward. For 
n = terms we have 



(0) 



16 



dr r 



2 2\R 



E^ 



and 



(0) 



16 



dr r 



Hp/R) 

\n{r/R) 



2 2 Vi? 



(B27) 



2 2\R 



2^+2R^- 2p2 _ [ ln(p/i?) - ^ (^)lj28) 



Using this and eq. (|B29p in eq. (|B34p we have, 



Eself /L 



327r 



n2 . 2 



(B34) 



Note that both the energy and force on the disclination 
due to the self-interaction vanish at the boundary as p — > 
R. 

Because a^^ contains only an monopole (n = 0) term, 
it is relatively straightforward to compute the mutual 
screening of the twist- and defect-induced stresses. 



E. 



(0) 



E 



i°)=8p^ln(p/i?) + 2i?^-2^. 



(B29) 



For n > 1 and E]^ ' we have. 



Es 



E. 



(n>l) 
< 



16 



dr r 



2n+l 



np" 





„ pn+2 







r dr r <j^'^\r)^{R^ -2r^)KT 
8 Jo 8 



(B35) 



(n + 1) 



1 (n + 1) / p \ 2™ / p \ 2nH 

n \r) ~ \r) 



(B30) 



For E^ 



(1) 



9 9 
r'^p^ 



i?4 



2- 



(DTI 



:8{2p2ln(i?/p)-2(i?2-p2)( 



{R'-P')(P 



P_ 

R 

2n 2 



2 r 



2i?2 



(Dhd)]} 



(n) 

and for n > 2 and i?> we have 



(ri>2) 



= 16 



r2"p2" [(n-l- 1) 



^4n 



1 ^p^2n 2 /p\2«r(n + 1) 
n \RJ L n 



. n'^in — 1) 



2{ P 



2n- 



--(i?^-p- 



Hi) 



R 

2"r(n+ 1) 



2n+2 _ p2n+2 \^ . ^ , 



1) 17?^* 



n 

2n 



{n + 1) 



Disclination interactions 



Here, we use the multipole expansion of the single 
disclination stress profile derived above, eqs. (|B19[) - 
(IB23[) . to calculate the elastic energy of interaction be- 
tween pairs of disclinations in the cylindrical cross section 
of the two-dimensionally order material. We consider two 
disclinations with respective topological charge, si and 
S2, located at polar coordinates (pi,0) and {p2,4>), and 
without loss of generality, assume p2 > pi. The elastic 
energy is computed from the following area integral. 



E,2/L = K,' J dA{4l + alDiatl + (B36) 



where ct^^ and cr^^ denotes the direct and induced 
stresses generated by the ath defect. Expanding the 
2 1 nroduct in the integrand above, we consider each con- 
tribution in turn, beginning with the induced-induced 
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terms. 

Eii — 



V (87r)2 J 

= 16^ ...{(lni?+i + J,)(lni?+i + J, 
^ ^rjpipzy + 1) _ ^pi j 



n=l 

(n + 1) 



cos 

,2 



M)} 



1 1 /PiP2\"r (n + l) /pi 



(n + 1) 



is(n(/))| 



+ 



-ln(l-^.^)+L.(^.^) 



i?2 
PlP2 



ln(l-^e^^ 



2 2 

i?4 



In 1 



PlP2 
i?2 ■ 



(B37) 



Next, we consider the direct-direct terms, 

'2„. „_ , _i /.fl 



" V (87r)2 ) i 

= 16|^'"rfr r[(lnpi + l)(lnp2 + l) 
+ -y^^( ) cos(n0) 

n=l 

+ J dr r[(lnr + l)(lnp2 + l) 
+ / dr r[(lnr + l)(lnr + l) 

n=l 

= 4|i?2 + p2 + (p2 ^ ^ 2i?2(l + Ini?) Ini? 



cosine 



■(-)' 



Ln2(n + 1) Vp2 

OO ^ 

, 1 r„2/PiA" T,i(PWiY 



^P2' 



+2pip2 ln(i?/p2) cos ( 



(B38) 



Making use of the series identities we have. 
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dd 



4Re 



i?2 + [pI +pl)\np2 + 2i?2 (1 + In i?) In i? - 2pip2 cos In (l - ^e'"^) + {pi + pl)\n (l 



pip^e^nn ( 1 - ^e'^ 



i?'Li2 



P2 



Pi 
— f 

P2 



i?2 



2pip2 ln(i?/p2)e'* (B39) 



r 



Finally, we calculate the induced-direct terms, 



Eid 







\ (87r)2 ) 


Jo 



dr r al'^{r)att{r) 



{^"d..[(lni?+l + ^)(lup2 + l 

^ OO 1 9 

1^ l/ Hpi X"/ 



2 ■ 2i?2y' 

„2 



cos(n(/i) 



P2 
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The corresponding term, Efj , is obtained by interchang- 
ing pi and p2 in Combining both induced-direct 
terms. 
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(B40) Evaluating the series, we have, 
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Combining eq. (|B37p . (|B39I) and (|B42I) we arrive at an 
expression for the elastic interaction between disclina- 
tions, Eint = LKos'^/{32Tr){Eu + Em + E,d) 
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}, (B43) 



Note that when either defect is pushed to the boundary, 
when pi — )• i? or p2 — J' R, |fi2 — = |ri — r2p and 
according to eq. (|B43p . Eint{j^i,J^2) = 0. Also, note that 
in the limit that 2 defects of charge si and S2 fuse and 

ri =r2, S„t(ri,ri)/L = i^oSiS2i?V(lM(l-P?/^')', 
which combined with the sum of the self-energies of the 
two elementary charges, yields the self energy of a discli- 
nation of total charge si + S2- 



where 



|fi2 - R|2 . i?2 ^ pM 



R^il 



R-^J 



2piP2 cos < 

d\ 



R^J 



I"! - r2 



|fB44) 



20 



[1] D. R. Nelson, Defects and Geometry in Condensed Matter [35] 
Physics, (Cambridge, Cambridge, 2002). 

A. A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). [36 
J.-F. Sadoc and R. Mosseri, Geometrical Frustration, [37 
(Cambridge, Cambridge, 1999). 

D. C. Wrigiit and N. D. Mermin, Rev. Mod. Piiys. 61, [38 
385 (1989). 

S. R. Renn and T. C. Lubensky, Phys. Rev. A 38, 2132 [39 
(1988). 

T. C. Lubensky and J. Prost, J. Phys. II France 2, 371 [40 
(1992). 

F. C. MacKintosh and T. C. Lubensky, Phys. Rev. Lett. 
V. Vitelh and A. M. Turner, Phys. Rev. Lett. 93, 215301 [41 
(2004). 

D. R. Nelson and L. Peliti, J. Phys. France 48, 1085 [42 

(1987) . [43 

C. D. Santangelo, V. Vitelli, R. D. Kamien and D. R. 
Nelson, Phys. Rev. Lett. 99, 017801 (2007). [44 

D. R. Nelson, R. D. Kamien, C. D. Santangelo and V. 
Vitelh, Phys. Rev. E 80, 051703 (2009). [45 
H. S. Seung and D. R. Nelson, Phys. Rev. A 38, 1005 

(1988) . [46 
M. J. Bowick, D. R. Nelson and A. Travesset, Phys. Rev. [47 
B 62, 8738 (2000). 

M. Bowick, A. Caccuito, D. R. Nelson and A. Travesset, [48 
Phys. Rev. Lett. 89, 185502 (2002). [49 
A. R. Bausch, M. J. Bowick, A. Cacciuto, A. D. Dins- 
more, M. F. Hsu, D. R. Nelson, M. G. Nikolaides, A. [50 
Travesset and D. A. Weitz, Science 299, 1716 (2003). 
W. T. Irvine, V. Vitelh and P. M. Chaikin, Nature 468, [51 
947 (2011). 

R. D. Kamien and D. R. Nelson, Phys. Rev. Lett. 74, [52 
2499 (1995). 

R. D. Kamien and D. R. Nelson, Phys. Rev. E 53, 650 [53 
(1996). 

T. J. Wess, in Gollagen: Structure and Mechanics, edited [54 
by P. Fratzl (Springer, Boston, 2008), Chap. 3. 
J. W. Weisel, J. Thromb. Haemost. 5 (Suppl. 1), 116 [55 
(2007). 

M. M. A. E. Claessens, C. Semmrich, L. Ramos and A. [56 
R. Bausch, Proc. Natl. Acad. Sci. USA 105, 8819 (2008). 
H. Shin, K. R. Purdy Drew, J. R. Bartles, G. C. L Wong [57 
and G. M. Grason, Phys. Rev. Lett 103, 138102 (2009). [58 
Y. Bouligand, C. R. Chimie 11, 281 (2008). 
J. W. Weisel, C. Nagawami and L. Makowski, Proc. Natl. [59 
Acad. Sci. USA 84, 8991 (1987). [60 
M. S. Turner, R. W. Briehl, F. A. Ferrone and R. Josephs, [61 
Phys. Rev. Lett. 90, 128103 (2003). [62 

G. M. Grason and R. F. Bruinsma, Phys. Rev. Lett. 99, 
098101 (2007). [63] 
G. M. Grason, Phys. Rev. E 79, 041919 (2009). 
Y Yang, R. B. Meyer and M. F. Hagan, Phys. Revl Lett. 
104, 258102 (2010). 

C. Heussinger and G. M. Grason 135, 035104 (2011). 
G. M. Grason , Phys. Rev. Lett. 105, 045502 (2010). 
A. Azadi and G. M. Grason, to be published (2011). 
O. StenuU and T. C. Lubensky, in preparation. 
J. V. Selinger and R. F. Bruinsma, Phys. Rev. A 43, 2910 
(1991). 

[34] P. G. de Gennes and J. Prost, The Physics of Liquid [64] 
Crystals (Claredon, Oxford, 1993), 2nd ed. 



G. Grinstein and R. A. Pelcovits, Phys. Rev. Lett. 95, 
856 (1981). 

M. Kleman and P. Ostwald, J. Phys. 47, 655 (1982). 
L. D. Landau and E. M. Lifshitz, Theory of Elasticity, 
3rd ed. (Pergamon, Oxford, 1986), Chap. 2. 
M.-F. Achard, M. Kleman, Y. A. Nastishin and H.-T. 
Nguyen, Eur. Phys. J. E 16, 37 (2005). 

A. A. Kornyshev, D. J. Lee, S. Leikin and A. Wynveen, 
Rev. Mod. Phys. 79, 943 (2007). 

P. M. Chaikin and T. C. Lubensky, Principles of Con- 
densed Matter Physics, Chap. 6 (Cambridge University, 
Cambridge, 1995). 

E. L. Starostin, J. Phys. Condens. Matter 18, S187 

(2006) . 

R. D. Kamien, Rev. Mod. Phys. 74, 953 (2002). 

H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. Curl and 
R. W. Smalley, Nature, 318, 162 (1985). 

D. L. D. Caspar and A. King, Cold Spring Harbor Symp. 
Quant. Biol. 27, 1 (1962). 

T. Aste and D. L. Weaire, The Pursuit of Perfect Pack- 
ing, 2nd ed. (Taylor and Francis, New York, 2008). 
P.-G. de Gennes, Solid State Commun. 10, 753 (1972). 
L. Giomi and M. Bowick, Phys. Rev. B 76, 054106 

(2007) . 

W. F. Brown, Phys. Rev. 60, 139 (1941). 

V. M. Bedanov and F. M. Peelers, Phys. Rev. B 49, 2667 

(1994). 

B. Partoens and P. Singha Deo, Phys. Rev. B 69, 245415 
(2004). 

M. Kong, B. Partoens, A. Matulis and F. M. Peelers, 
Phys. Rev. E 69, 036412 (2004). 

R. Bubeck, C. Bechinger, S. Neser and P. Leiderer, Phys. 
Rev. Lett. 82, 3364 (1999). 

A. V. Filinov, M. Bonitz and Y. E. Lozovik, Phys. Rev. 
Lett. 86, 3851 (2001). 

B. J. Baelus, F. M. Peelers and V. A. Schweigert, Phys. 
Rev. B 63, 144517 (2001). 

H. Saarikoski, A. Harju, M. J. Puska and R. M. Niemi- 
nen, Phys. Rev. Lett. 93, 116802 (2004). 
See e.g., L. Bonsall and A. A. Maradudin, Phys. Rev. B 
15, 1959 (1977). 

R. deWit, J. Appl. Phys. 42, 3304 (1971). 

W. F. Harris and L. E. Scriven, J. Appl. Phys. 42, 3309 

(1971). 

R. D. Kamien, J. Phys. II France 6, 461 (1996). 

M. Kleman, J. Phys. Lett. 46, L-723 (1985). 

M. Kleman, J. Phys. 46, 1193 (1985). 

T. Yue and X. Zhang, J. Phys. Chem. B. 115, 11566 

(2011). 

The ring of 5-fold defects in Fig. [IJa is reminiscent of the 
disc-packing textures on spherical surfaces studied in M. 
Rubenstein and D. R. Nelson, Phys. Rev. B 28, 6377 
(1983). Such packings were constructed by determinis- 
tic algorithm that grows locally-dense clusters radially 
outward from a given pole of the sphere. While regions 
small compared to sphere radius retain hexagonal pack- 
ing, upon reaching a critical lateral radius, these pack- 
ings exhibit characteristic regions rich with excess 5-fold 
disclinations. 

It is intriguing to note that both the theory of double- 
twisted A'' -|- 6 order [H^ and the model of ideal, twisted 
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filament packing in [goI . [gH find that both problems 
map on to the ordering on a sphere of curvature Kg ~ 
Sl^ , that is 3 times smaller than the intrinsic twist result 



of eq. (|2ip for twisted and 2D order filament packings in 
Euclidean space. 



